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Abstract 

Background: Lawsonia intracellularis is an obligate intracellular bacteriunn and the etiologic agent of proliferative 
enteropathy. The disease is endemic in pigs, emerging in horses and has been described in various other species 
including nonhuman primates. Cell proliferation is associated with bacterial replication in enterocyte cytoplasm, but 
the molecular basis of the host-pathogen interaction is unknown. We used laser capture microdissection coupled 
with RNA-seq technology to characterize the transcriptional responses of infected enterocytes and the host- 
pathogen interaction. 

Results: Proliferative enterocytes was associated with activation of transcription, protein biosynthesis and genes 
acting on the Gi phase of the host cell cycle (Rho family). The lack of differentiation in infected enterocytes was 
demonstrated by the repression of membrane transporters related to nutrient acquisition. The activation of the 
copper uptake transporter by infected enterocytes was associated with high expression of the Zn/Cu superoxide 
dismutase by L intracellularis. This suggests that the intracellular bacteria incorporate intracytoplasmic copper and 
express a sophisticated mechanism to cope with oxidative stress. 

Conclusions: The feasibility of coupling microdissection and RNA-seq was demonstrated by characterizing the 
host-bacterial interactions from a specific cell type in a heterogeneous tissue. High expression of L intracellularis 
genes encoding hypothetical proteins and activation of host Rho genes infers the role of unrecognized bacterial 
cyclomodulins in the pathogenesis of proliferative enteropathy. 
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has been associated with carcinogenesis in chronic dis- rodentium) are able to promote enterocyte proliferation 

eases caused by Helicobacter pylori, Salmonella typhi with minimal inflammatory responses [4,5]. C. roden- 

and Citrobacter rodentium [1]. Proliferative changes tium induces proliferation of mouse colonic enterocytes 

resulting in a hyperplastic but non-carcinogenic process when it attaches to the apical membrane [6]. It is be- 

is induced in Lawsonia intracellularis and Bartonella lieved that L intracellularis escapes from the endosome 

spp. infections [1]. The inflammatory mediators gener- after internalization and multiplies freely in the cyto- 

ated during the chronic gastritis in H, pylori infections plasm of undifferentiated intestinal crypts promoting 

have been related with oxidative DNA damage and cell their proliferation and progressive replacement of the 

differentiated intestinal epithelium by immature infected 
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infection in crypt cells [7]. These unusual pathological 
changes characterized by the presence of a large number 
of intracellular bacteria and proliferation of enterocytes 
suggests that L, intracellular is has adopted mechanisms 
of survival and pathogenesis that are unique among bac- 
terial pathogens [5]. To date, hypotheses and specula- 
tions have been discussed regarding the pathogenesis of 
this infection but, the underlying mechanisms by which 
L, intracellularis induces proliferative changes have not 
been addressed. 

Lawsonia intracellularis is an obligate intracellular 
bacterium and the etiologic agent of proliferative enter- 
opathy (PE) [8]. Mild to severe diarrhea is the major 
clinical sign described in infected animals and is directly 
associated with cell proliferation and replication of the 
bacteria in the intestinal epithelium [9]. PE is endemic in 
swine herds, an emerging disease in horses and has been 
reported in various other species, including non-human 
primates, wild mammals and ratite birds [10-13]. Micro- 
scopically, different levels of hyperplasia in the intestinal 
crypts associated with significant reduction in the num- 
ber of goblet cells are observed in animals affected with 
PE. Although inflammation is often observed at late 
stages of the disease, it does not represent a primary 
characteristic associated with L, intracellularis infection. 

The disease was first reported in 1931 but, due to its 
fastidious properties, porcine L, intracellularis was first 
isolated in 1993 using rat small intestinal cells cultu- 
red in a strictly defined microaerophilic environment 
[14,15]. Since then, the chronological dynamics of infec- 
tion in vitro and in vivo have been well-characterized 
[16,17]. While proliferative changes in vivo follow an in- 
crease in the number of intracellular bacteria, the bacter- 
ium is not able to induce proliferation in infected cells 
in vitro [18]. This intriguing observation associated with 
the fastidious properties of this bacterium has limited 
the studies on the pathogenesis of L, intracellularis. 
Additionally, the adaptation of this microaerophilic, but 
obligate intracellular organism to grow freely in the 
cytoplasm of metabolically active enterocytes suggests 
that this organism has properties that are novel and 
unique among bacterial pathogens. 

In the present study, we hypothesized that genes 
expressed by L intracellularis in the host cytoplasm are 
capable of inducing proliferation and preventing diffe- 
rentiation of immature enterocytes by altering cell 
cycle-associated pathways. We established a method in- 
tegrating laser capture microdissection (LCM) and 
RNA-seq technology to characterize the host gene ex- 
pression and the host-bacterial interactions in vivo. Acti- 
vation of transcription, protein synthesis and Rho family 
genes in infected enterocytes characterized the transcrip- 
tional mechanisms involved in the cell proliferation. The 
ability of L, intracellularis in preventing enterocyte 



differentiation and maturation was proved by the con- 
sistent down-regulation of apical membrane transporters 
related to nutrient acquisition by infected enterocytes. 
The intracellular bacteria showed high level of expres- 
sion of a sophisticated oxidative protection mechanism 
which involves redox enzymes and a rubrerythrin- 
rubredoxin operon (rubY-rubA), Rho genes expressed by 
the host and associated with the high expression of bac- 
terial genes encoding hypothetical proteins implies a 
potential role of unrecognized bacterial effector proteins 
that modulate eukaryotic cell cycle (cyclomodulins) in 
the pathogenesis of PE. 

Results 

Technical characterization of the LCM coupled with RNA- 
seq analysis 

Based on previous studies reporting the chronological 
course of the L intracellularis infection in vivo [19,20], 
experimentally-infected pigs in our study were moni- 
tored every other day regarding clinical signs (diarrhea 
and general attitude) and quantitative fecal shedding of 
L intracellularis DNA until the end point of the study 
(21 PI). The negative control group was also monitored 
and its negative status was consistent throughout the 
study period. One animal from the infected group 
exhibiting only mild to moderate lesions typical of PE 
associated with grade 1 and 2 in the IHC grading 
scheme was discarded from the experiment. In order to 
have the same number of replicates in both infected and 
control groups, one animal was randomly discarded 
from the negative control group. As a result, five bio- 
logical replicates (animals) from each group were in- 
cluded in the LCM procedures and RNA-seq analysis. 

The quality of RNA samples extracted from microdis- 
sected tissues was evaluated by assessing bacterial and 
eukaryotic ribosomal RNAs using Agilent Bioanalyzer 
2100 (Figure IC). The successful recovery of high quality 
eukaryotic RNA from microdissected cells using the LCM 
procedures described in our study has been well- 
established [21]. Since studies using prokaryotic RNA 
from microdissected cells have not been well-explored, the 
feasibility of using L intracellularis mRNA for sequencing 
was confirmed by one-step RT-PCR using specific primers 
targeting three housekeeping genes (Figure ID). 

The sequence reads representing the RNA transcripts 
derived from the host cells and the intracellular bacteria 
were mapped onto both pig {S. scrofa 10.2) and the 
L intracellularis (PHE/MNl-00) reference genomes 
available at NCBI. From a total average of 22,136,064 
reads generated, 82% (18,280,138) and 4% (1,000,440) 
were mapped against the porcine and the bacterium 
genome, respectively (Additional file 1: Table S2). A total 
of 11,778 genes were expressed in the porcine genome 
and met the criteria for differential expression analysis 
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Figure 1 Laser capture microdissection of an intestinal crypt infected with L Intracellularis (A-B) and evaluation of RNA quality from 
microdissected cells (C-D). (A) Hematoxylin stained cryosection of infected ileal mucosa. (B) Microdissected intestinal crypt captured in the 
thermoplastic film of the LCM cap prior to RNA isolation. (C) Agilent Bioanalyzer data showing bacterial and eukaryotic ribosomal RNA in infected 
and non-infected cells. (D) One-step RT-PCR products of three protein-encoding genes of L. intracellularis. 



(see Materials and Methods). Up-regulation of 119 and 
down-regulation of 46 protein coding genes were identi- 
fied in infected enterocytes (fold-change >2.0; p<0,05) 
(Additional file 2: Table S3). 

In the bacterial gene expression analysis, 754 protein 
coding genes had at least one mapped read against the 
reference DNA sequence of the L intracellularis PHE/ 
MNl-00 isolate, which was the same isolate used in 
the experimental infection. The top 20 genes highest 
expressed by L, intracellularis and their respective 



normalized transcription levels based on the RPKM 
values are described in Table 1. 

In order to screen for the presence of any other micro- 
bial organisms in the experimental samples, the se- 
quence reads were also mapped against viral, bacterial 
and archaeal genome databases. While no detectable hits 
were found against the viral and archaeal genomes, 
0.052% and 0.084% of the sequence reads from the 
infected and control group, respectively, were mapped 
against the bacterial database. Three unculturable 
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Table 1 The top 20 highly expressed genes by L intracellularis In the enterocyte cytoplasm 



Locus 


Gene product 


Predicted function 


Log2 (RPKM) 


LI0461 


Hypothetical protein 


Unknown 


15.0 


LI0625 


Chaperonin GroEL 


Protein folding 


14.9 


Lll 159 


Hypothetical protein 


Unknown 


13.7 


LI0447 


Hypothetical protein 


Unknown 


13.6 


LI0935 


Elongation factor Tu {tufA) 


Protein biosynthesis 


13.5 


LI0624 


Co-chaperonin GroES 


Protein folding 


13.5 


LI0005 


Cu-Zn superoxide dismutase precursor {sodQ 


Oxidative stress 


13.3 


LIC060 


Hypothetical protein 


Unknown 


12.9 


LI0267 


Hypothetical protein 


Unknown 


12.8 


LI 1075 


Dioxygenases related to 2-nitropropane dioxygenase 


Oxidative stress 


12.7 


LI0809 


Hypothetical protein 


Unknown 


12.5 


LI0912 


Molecular chaperone DnaK 


Protein folding 


12.4 


LI0439 


Hydrogenase-1 small subunit 


Cell metabolism 


12.4 


LI0902 


Outer membrane protein related to OmpA-OmpF porin 


Membrane transport 


12.4 


Lll 158 


Hypothetical protein related to secretion system effector C 


Intracellular survival 


12.3 


LI 1005 


Pseudouridine synthase [truB) 


Protein biosynthesis 


12.2 


LI0559 


50S ribosomal protein LI 3 


Protein biosynthesis 


12.2 


LI0697 


Rubrerythrin {rubY) 


Oxidative stress 


12.2 


LI0043 


Hypothetical protein 


Unknown 


12.0 


LIC103 


Methyl-accepting chemotaxis protein {pilJ) 


Chemotaxis 


11.9 



bacterial species {Candidatus Zinderia, Candidatus 
Carsonella ruddii, Candidatus Sulcia muelleri) and 
Propionibacterium acnes were commonly found in 
both groups. An uncharacterized species of the genus 
Peptoniphilus previously reported in the oral cavity of a 
human was exclusively found in the control group. 

The RNA-seq expression data were validated by qRT- 
PGR based on the relative quantification of 16 differen- 
tially expressed genes from the host and 10 bacterial 
genes (Additional file 3: Table SI). The averages of fold- 
change in gene expression from infected enterocytes and 
the relative transcriptional levels of L, intracellularis 
genes were plotted against the log2 (RPKM) (Figure 2). 
The log2 transformed fold-change from the host gene 
and transcript levels from the bacterial genes were posi- 
tively correlated on a linear regression model (/7-value < 
0.05). Based on the qRT-PCR validation, the RNA-seq 
data properly estimated the fold-change expression in 
infected enterocytes and the transcription levels of L, 
intracellularis genes. 

Protein biosynthesis and transcription 

The 165 porcine genes differentially expressed in Lawso- 
nia-infected enterocytes were analyzed regarding their 
biological functions and molecular networks based on 
the mammalian gene expression information available in 
the Ingenuity® System. The system knowledge database 
recognized and analyzed 144 differentially expressed 
genes. The IPA system associates the set of differentially 



expressed genes with cellular networks (focus genes) and 
creates a score based on the number of network eligible 
genes they contain (Table 2). The protein biosynthesis 
network was most correlated with the set of genes differ- 
entially expressed and up-regulated in LawsoniaAnfected 
enterocytes. Figure 3A shows the molecular interaction 
of this network which mainly includes ribosomal pro- 
teins and mRNA translation factors. The eukaryotic 
initiation factor (IEF2) signaling was the canonical cell 
pathway (CP) more significantly associated (j^-value 
2.76E-22) with the genes differentially expressed 
(Figure 3B). A remarkable influence of IEF2 signaling on 
the protein biosynthesis was identified by merging both 
network and pathway analysis (Figure 3A). 

Functional clustering analysis using DAVID knowledge 
database confirmed the involvement of the 27 differen- 
tially expressed genes in the protein biosynthesis (data not 
shown) and also showed the association of four genes act- 
ing in the positive regulation of transcription (Table 3). 
The activation of these two cellular processes (transcrip- 
tion and protein biosynthesis) revealed a global increasing 
in the cell metabolism in response to the L intracellularis 
infection and it has also been described to occur during 
the gap phase 1 (Gi) of the host cell cycle [22]. 

Cell cycle and apoptosis 

In addition to the positive regulation of transcription, 
DAVID clustering analysis identified genes associated 
with cell cycle and apoptotic events differentially 
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Figure 2 Correlation between RNA-seq and qRT-PCR expression data. (A) Plot of the relative fold-change in gene expression of 16 porcine 
genes from infected enterocytes by RNA-seq (x-axis) and qRT-PCR (y-axis). R^ = 0.92 (p<0.05). (B) Plot of the relative quantification often L 
introcelluloris genes by RNA-seq (x-axis) and the transcript levels generated by qRT-PCR (y-axis). R^ = 0.81 (p<0.05). 



expressed in infected enterocytes (Table 3). Cell 
cycle-associated genes mainly represented by Ras 
homolog proteins were significantly induced in associ- 
ation with Cyclin-dependent kinase 2 (CDK2), The 
aberrant activation of Rho-genes including those 



described in our study has been well-described in 
oncogenesis by causing deregulation of cell cycle 
progression and promoting cell proliferation [23]. 
Additionally, the co-expression of Rho proteins and 
Cyclin-dependent kinases (CDKs) specifically act by 



Table 2 The four cell networks most associated with the genes differentially expressed in /.aui^son/a-infected 
enterocytes 



Associated network functions 


Score* 


Focus genes 


Protein Synthesis, Cellular Assembly and Organization, Molecular Transport 


59 


27 


RNA Post-Transcriptional Modification, Cancer, Hematological Disease 


42 


21 


Cellular Movement, Cancer, Reproductive System Disease 


31 


17 


Cell Death and Survival, Carbohydrate Metabolism, Molecular Transport 


31 


18 


* Based on the number of network eligible genes within the referred network. 
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Figure 3 Cellular network and canonical pathway analysis of genes differentially expressed by enterocytes infected with L 
intracellularis. (A) Molecular interaction representing part of the protein biosynthesis network which was most correlated with the set of genes 
significantly up-regulated in infected cells. Strong interaction between the main canonical pathway (CP: EIF2 signaling) identified in the 
differentially expressed genes and the protein biosynthesis network. (B) Canonical pathways (x-axis) most associated with genes differentially 
expressed based on the -log of p-value calculated by Fisher's exact test (y-axis). 
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Table 3 Functional clustering associated with genes differentially expressed in Laui^son/a-i nfected enterocytes 



Gene Symbol 


Gene product 


ln(fold 


Clustering 




change) 


Positive regulation of transcription 


TEFl 


Transcriptional enliancer factor 1 


2.2 


SMARCC2 


SWI/SNF complex subunit SMARCC2-like 


2.3 


SOX-9 


Transcription factor SOX-9 


2.6 


FHL2 


Four and a half LIM domains 2 


2.1 


Mitotic cell cycle-associated genes 


CDK2 


Cyclin-dependent kinase 2 


2.6 


RHOA 


Ras homolog family member A 


2.2 


RHOB 


Ras homolog family member B 


2.2 


RRAGA 


Ras-related GTP-binding protein A 


2.3 


RPL24 


60S ribosomal protein L24-like 


2.7 


SKA2 


Spindle and kinetochore-associated protein 2-like 


2.1 


PSMA1 


Proteasome subunit alpha type-l-like 


2.4 


Pro-apoptosis-related genes 


SLAM7 


Signaling lymphocytic activation molecule 7 


-4.6 


C6 


Complement component C6 


2.3 


C9 


Complement component C9 


3.5 


PRDXl 


Peroxiredoxin 1 


2.2 


SST 


Somatostatin 


-3.4 


TNFSFIO 


Tumor necrosis factor ligand 10 


-2.5 


VAV2 


Vav 2 guanine nucleotide exchange factor 


-2.1 


Anti-apoptosis-related genes 


DADl 


Defender against cell death 1 


2.3 


NME2 


Non-metastatic cells protein 


2.9 


PRNP 


Prion protein 


2.6 


PPP2CB 


Protein phosphatase 2, catalytic subunit, beta isozyme 


2.1 


RHOA 


Ras homolog family member A 


2.2 


TPTl 


Translationally-controlled tumor protein-like 


2.7 


ERBB3 


Receptor erythroblastic leukemia viral oncogene 


2.2 



stimulating the entry and progression of the Gi phase 
of the host cell cycle [24]. 

The functional clustering analysis revealed 14 differ- 
entially expressed genes associated with pro- and 
anti-apoptotic events (Table 3). All seven genes re- 
lated to anti-apoptotic events were up-regulated in 
infected cells. Among the pro-apoptotic genes, four 
were down-regulated and three up-regulated. Interest- 
ingly, two of these pro-apoptotic genes that are also 
involved in the cellular immune response against 
intracellular pathogens (Signaling lymphocytic activa- 
tion molecule 7 and Tumor necrosis factor ligand 10) 
were down-regulated. On the other hand, significant 
up-regulation of genes encoding the major histocom- 
patibility complex class I (MHC-I) was identified in 
LawsoniaAnfected enterocytes (Figure 4A). 



Nutrient acquisition and electrolyte secretion 

Consistent down-regulation of numerous genes expres- 
sed in the apical membrane of enterocytes that are 
involved in nutrient acquisition was observed in 
Lawsonia-infected cells (Figure 4A). These membrane 
transporters are involved in the absorption of carbohy- 
drates (sodium/glucose co-transporter and sucrase- 
isomaltase), amino acids (cationic amino acid trans- 
porter), bile acid (sodium/bile acid co-transporter), lipids 
(lipid phosphate phosphohydrolase) and Vitamin B12 
(cubilin receptor). Additionally, the intracellular infec- 
tion also affects the electrolyte secretion by decreasing 
expression of the chloride channel gene (CLCAl). The 
reduction of both nutrient acquisition and electrolyte se- 
cretion indicates that L intracellularis may be able to 
prevent cell differentiation in immature enterocytes. 
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Figure 4 Proposed model for host-pathogen interaction in porcine enterocytes infected with L intracellularis. (A) Infected enterocyte. 
Apical membrane exhibiting down-regulation of genes involved in nutrient acquisition and electrolyte secretion (green) and up-regulation of 
copper uptake protein (red). Basolateral membrane exhibiting up-regulation of glucose transporter (SLC2A1) and MHC class I genes. (B) 
Intracellular bacterium. High expression (yellow) of genes included in the oxidative stress protection system: redox enzymes, Cu-Zn superoxide 
dismutase {sodQ and rubrerythrin-rubredoxin {rubY-rubA) operon (Adapted from Lumppio et al. 2001). Moderate expression genes encoding the 
basal body of the type III secretion system (TIIISS) (light yellow) and high expression of TIIISS effector proteins {PcrH-SseC-Sse? operon). Genetic 
organization and gene expression (log2 [rpkm]) of the rubY-rubA and PcrH-SseC-Sse operons. 



In contrast to the down-regulation of genes expressed 
on the apical membrane, the gene encoding the glucose 
transporter 1 (also known as solute carrier family 2, glu- 
cose transporter member 1) was highly up-regulated in 
infected enterocytes. This transporter was the third 
highest expressed gene (Additional file 2: Table S3) and 
its expression also has been reported on the basolateral 
membrane of human enterocytes in vitro [25] and rat 
jejunum in vivo [26]. Additionally, significant up- 
regulation of the high-affinity copper uptake protein 
(CTRl) involved in copper absorption was also found in 
infected enterocytes (Figure 4A). 

Host-pathogen interaction 

The transcriptional landscape of the intracellular bac- 
teria was determined by classifying all genes with at least 
one mapped read into one of three levels of expression: 
low, moderate and high. As expected, genes encoding 
ribosomal-related proteins were the functional category 
most associated with the mapped reads and the majority 
of these genes exhibited moderate to high expression. 
The L intracellularis gene expression was characterized 
and discussed based on the most highly expressed bac- 
terial genes (Table 1). 



The functional categories of bacterial genes that were 
highly expressed included those involved in protein fold- 
ing (e.g. groES'groEL operon and chaperone dnaK) and 
biosynthesis (e.g. elongation factor Tu), oxidative stress 
(e.g. Cu-Zn superoxide dismutase, rubY-rubA operon 
and dioxygenases), secretion system effector- related pro- 
teins {PcrH'SseC'llllS9 operon) and various hypothet- 
ical proteins. Since we were not able to distinguish free 
bacteria in the cell cytoplasm from those organisms 
within cell endosomes, the host-pathogen interaction was 
discussed considering both scenarios. Figure 4B illustrates 
the biological activities of the bacterial proteins encoded 
by highly expressed genes considering the intracellular 
microenvironment of Lawsonia-infected enterocytes. 

While redox enzymes catalyze reduction of O2 derived 
from the cell cytoplasm, Cu-Zn superoxide dismutase C 
(sodC) and rubrerythrin-rubredoxin operon (rubY-rubA) 
neutralize reactive oxygen species (O2 and H2O2) gener- 
ated in the endosome or from the reduced O2 molecule 
(Figure 4B) [27]. An operon genetically related to the 
Salmonella pathogenicity island 2 (SPI2) composed of 
three genes, a chaperone {PcrH), an effector protein {SseC) 
and a hypothetical protein LI1159 (referred as Sse?) was 
highly expressed (Figure 4B). Interestingly, genes encoding 
the type III secretion system (TIIISS) apparatus whereby 
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these effector proteins would be delivered to the cell cyto- 
plasm showed only moderate expression (data not shown) 
and were located downstream of the SPI2-related operon 
on the L. intracellularis chromosome. 

Among the ten genes most highly expressed by the 
intracellular bacterium, five represented hypothetical 
proteins. A summary of the predictive analyses evaluat- 
ing their structures and biological functions are de- 
scribed in the supplemental material (Additional file 4: 
Table S4). Two different families of proteins were pre- 
dicted for the gene locus LI0447 (aminomethyltrans- 
ferase beta-barrel and growth factor receptor domain). 
Additionally, transmembrane proteins (porin and 
autotransporter), proline isomerase, extracellular prote- 
ase and secretory factor were also identified based on 
predictive motifs. 

Discussion 

The present study used LCM technology to microscop- 
ically dissect enterocytes from pigs experimentally 
infected with L intracellularis and to characterize the 
cell and the bacterial gene expressions using high 
throughput sequencing (RNA-seq). The results demon- 
strated the usefulness of coupling LCM and RNA-seq 
techniques to study the host-pathogen interaction in a 
specific cell population present in a heterogeneous tis- 
sue. The simultaneous evaluation of the gene expression 
changes in both the host and the pathogen has been re- 
cently designated dual RNA-seq [28]. The challenge in 
studying both sets of transcriptional profiles from a 
common sample arises mainly because of the extensive 
variety of the bacterial organisms regarding their gen- 
omic compositions (e.g. %CG contents) and the particu- 
lar characteristics of each infectious process at the 
cellular level. Since the abundance of eukaryotic RNA is 
significantly higher compared to the prokaryotic RNA in 
an infected cell, a sufficient number of bacteria per host 
cell followed by an unbiased RNA amplification step are 
crucial for studying the cell-bacteria interactions [28]. In 
light of this, we established the endpoint of our study 
(21 days PI), based on the chronological course of the L, 
intracellularis infection previously described [17,29]. 
This previous study demonstrated a progressive increase 
in the number of L. intracellularis in the feces of 
infected pigs from day two (10^) through day 22 (10^) 
post-infection. Confirming these previous observations, 
our study showed at 21 days PI, numerous bacteria in 
proliferative enterocytes indicating the active stage of 
the infectious process. We hypothesized that the disease 
at this point was approaching its peak when the intracel- 
lular bacteria were exhibiting their virulence factors and 
the host cells were responding to the pathogen at appro- 
priate levels to allow us to evaluate the changes in gene 
expression. We believe that earlier endpoints would not 



provide a sufficient amount of RNA to be recovered and 
later stages of the disease would not represent the loga- 
rithmic phase of the bacterial growth. Additionally, the 
increasing presence of L. intracellularis within the lam- 
ina propria along the disease progression would also re- 
duce the number of bacteria in enterocytes. In addition 
to the limitation of having only one time point, our 
study also could not evaluate the role of soluble proteins 
and growth factors secreted into the extracellular matrix 
of the lamina propria and their influence in the entero- 
cyte proliferation and differentiation. 

A previous chronological analysis using microarray 
technology evaluated the host response of fibroblastic 
cells in vitro at three time points (24, 48 and 72 hours 
PI) [30]. Altered transcription of genes related to cell 
cycle and cell differentiation were described [30]. How- 
ever, cellular proliferation which is the main phenotypic 
characteristic of the L, intracellularis infections in vivo 
has not been reproduced in in vitro models [14,18]. Fur- 
thermore, the transcriptional response of mesenchymal 
cells which are not the natural target cells for L, 
intracellularis needs to be interpreted with caution. 

A single time point analysis of the host gene expres- 
sion profile using intestinal tissues from pigs naturally 
affected with PE was recently described [31]. Although 
this study provided an interesting snapshot of the tran- 
scriptional host response, it used field cases of diarrhea 
where the samples were co-infected with porcine 
circovirus type 2 which may be a confounding factor in 
the evaluation of the expression of genes especially re- 
lated to the immune response [31]. Additionally, the 
microarray was performed using entire intestinal tissues. 
As a result, the specific characterization of the transcrip- 
tional host response was impaired by either the hetero- 
geneity of the cell population included in the intestinal 
tissues or the virus infection. By applying the LCM tech- 
nique to isolate ileal enterocytes in our study, we con- 
firmed the identification of several genes specifically 
expressed in the intestinal epithelium of the ileum (e.g. 
ileal sodium/bile acid cotransporter) corroborating with 
a recent comprehensive study describing the gene ex- 
pression atlas of the domestic pig [32]. 

Gene network and pathway analyses revealed that pro- 
tein biosynthesis and activation of transcription were the 
host cellular events most associated with the genes dif- 
ferentially expressed in Lawsonia-iniected enterocytes 
(Figure 3A). In agreement with these molecular findings, 
ultrastructural studies in vivo using electron microscopy 
identified the numerous bacteria occupying the apical 
cytoplasm of infected enterocytes which was otherwise 
composed almost entirely of free ribosomes and 
scattered mitochondria [29,33]. These morphological 
findings associated with our molecular results suggest a 
global increase in the cell metabolism in response to L 
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intracellularis infection. Taking the host cell cycle events 
into consideration, an increase in protein synthesis and 
transcription is required during the Gi phase in order to 
prepare the cell for subsequent division [22]. Interest- 
ingly, this protein synthesis network was also associated 
with cell proliferation at an early stage of infection with 
human immunodeficiency virus [34]. 

The differential expression of cell cycle-associated 
genes identified by the functional clustering analysis 
(Table 3) showed specific activation of Rho family genes 
(RhoA, RhoB and Rho GTPase). These molecules play a 
key role in carcinogenesis through their aberrant activa- 
tion which results in cell proliferation [23]. Rho family 
proteins specifically act on the Gi -checkpoint of the cell 
cycle when the transition to commit the cell to the pro- 
liferative stage occurs. If the signals responsible for pro- 
moting this transition are not present then the cells 
enter into the non-proliferative phase (Gq) [35]. Add- 
itionally, a gene encoding Rho GTPase was also highly 
up-regulated during L, intracellularis infection in vitro 
[30], suggesting that exacerbated activation of the Gi 
phase is an important mechanism involved in the prolif- 
erative changes induced by L, intracellularis in infected 
enterocytes. 

In addition to their roles in cancer development, Rho 
proteins are also pathologically activated by bacterial 
toxins also known as cyclomodulins [1,35,36]. Cytotoxic 
necrotizing factor (CNF) found in uropathogenic 
Escherichia coli, Pasteurella multocida toxin (PMT) and 
dermonecrotic toxin of Bordetella spp. act directly on 
Rho family proteins to bring about their irreversible acti- 
vation [37]. Furthermore, CNF-expressing E, coli estab- 
lish a persistent intracellular infection in the urogenital 
tract and suppress apoptosis by affecting the transcrip- 
tion levels of Bcl-2 family genes [37]. It is thought that 
apoptosis inhibition in target cells may favor bacterial 
persistence at the epithelium surface, thereby favoring 
bacterial replication and spread inside the host cell [38]. 
Although our study showed a predominant activation 
of anti-apoptotic-related genes compared with pro- 
apoptotic events, the Bcl-2 gene was not significantly 
activated in infected cells (log2-fold change = 1.06). 
Therefore, other mechanisms may be involved in the 
predominant activation of anti-apoptotic genes identified 
in the present study. 

The contribution of apoptotic mechanisms for the 
pathogenesis of L, intracellularis infections has been 
speculated over the years and still needs to be eluci- 
dated. Initially, a temporary reduction in apoptosis was 
hypothesized to be an important mechanism involved in 
the cell proliferation [29]. Later studies suggested an in- 
crease in apoptosis based on the Caspase-3 immunohis- 
tochemical staining [8,39]. One of these studies 
described the dynamic of the Caspase-3 staining through 



the chronological evaluation of experimentally-infected 
animals and showed variations on the pattern of 
Caspase-3 staining within different parts of the intestinal 
mucosa on day 19 PI [39]. Our study showed activation 
of the gene-encoding Caspase-3 (log2-fold change = 
1.12) in infected enterocytes, but the level did not reach 
the parameters for statistical significance. Considering 
the complexity and dynamism of the apoptosis process, 
the stage of the infection at the cellular level may dir- 
ectly influence the apoptotic gene network. Therefore, 
an in vitro model displaying the proliferative phenotype 
would be an ideal starting point to address this question. 
Despite the reduced expression of these two immune 
response-associated genes identified in our study, we ob- 
served significant up-regulation of MHC-I genes in 
infected cells, indicating that Lawsonia-dQviYQd antigen 
is presented to the lamina propria through the basal 
membrane of infected enterocytes (Figure 4A). 

Although the physiopathology of diarrhea in L, 
intracellularis infections remains to be elucidated [40], 
the significant down-regulation of numerous genes re- 
lated to nutrient acquisition observed in the present 
study indicates that malabsorptive diarrhea represents 
the major mechanism involved in the poor performance 
and growth of affected animals. Supporting our observa- 
tions, a lower absorption of glucose and electrolytes was 
reported using a hamster experimental model of PE [41]. 
The reduced expression of nutrient acquisition-related 
genes also indicates the lack of cell differentiation in 
infected enterocytes. 

The deficiency in nutrient acquisition described above 
seems to contrast with the increase in the cell metabol- 
ism characterized by the activation of transcription and 
protein synthesis-related genes also described here. As 
proposed in the Figure 4A, the high expression of the 
glucose transporter SLC2A1 on the basolateral mem- 
brane appears to compensate for this lack of nutrient ac- 
quisition from the intestinal lumen, since it has been 
described as being expressed in the basolateral mem- 
brane of human enterocytes in vitro [25] and rat je- 
junum in vivo [26]. 

The only gene significantly up-regulated in infected 
cells which is involved in nutrient acquisition and is 
physiologically expressed on the apical membrane is the 
high-affinity copper uptake gene [42]. Copper is an es- 
sential metal used by eukaryotic cells as a biochemical 
cofactor especially during the process of oxygen reduc- 
tion by cytochrome c oxidase which leads to the produc- 
tion of ATP [43]. The production of ATP by the 
eukaryotic host is crucial for metabolism of Lawsonia 
and two other groups of obligate intracellular bacteria 
(Chlamydiales and Rickettsiales), These organisms ex- 
press ATP/ADP translocase which catalyze the exchange 
of bacterial ADP for host ATP allowing bacteria to 
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exploit their hosts' energy pool, a process referred to as 
energy parasitism [44,45]. In addition to this essentiality 
for host cells, an increase in copper uptake has also been 
reported as a defense mechanism against intracellular 
pathogens because of its toxic properties. Free intracel- 
lular copper can lead to oxidative stress whereby cycling 
of copper oxidation and reduction produce reactive oxy- 
gen species through the Fenton reaction [43]. Interest- 
ingly, our study identified L, intracellularis genes related 
to oxidative stress highly expressed within the host cyto- 
plasm (Figure 4B). 

Three to five percent of the sequence reads generated 
in our study also mapped against the L, intracellularis 
genome. Less than 0.1% of the sequence reads mapped 
against four bacterial species commonly found in both 
infected and non-infected groups, when screened for the 
presence of virus and other bacteria. While three of 
these species represent unculturable organisms {Candi- 
datus) related to insect symbionts, one {Propionibacter- 
ium acnes) is part of the normal flora of the skin, oral 
cavity, large intestine, the conjunctiva and the external 
ear canal of humans [46]. The identification of the genus 
Peptoniphilus, found exclusively in the control group, 
may be due to human contamination during the experi- 
mental procedures, since it was reported in the oral 
cavity through the human microbiome project (NCBI 
accession: PRJNA52051). 

L, intracellularis genes belonging to the oxidative 
stress protection were the functional cluster most associ- 
ated with genes highly expressed in the host cytoplasm. 
Although this fastidious organism requires a strict 
microaerophilic environment for cultivation in cell cul- 
ture, Lawsonia is often located close to the cell mito- 
chondria, where the transport of oxygen is continuous 
due to the oxidative phosphorylation [47]. Additionally, 
this intracellular location is important to exploit the host 
energy pool by exchanging bacterial ADP for host gener- 
ated ATP, as discussed earlier. Based on this scenario, 
our study demonstrated that L intracellularis displays a 
sophisticated mechanism to survive in this microenvir- 
onment by coping with the oxidative stress and poten- 
tially incorporating intracytoplasmic copper taken by the 
host to combat the intracellular infection (Figure 4A). 

Our study identified high expression of an operon 
(PcrH'SseC'Sse) genetically related to the SPI2. The 
proteins included within this operon act as a translocon 
attached to the phagosomal membrane to allow trans- 
location of effector proteins into the cell cytoplasm [48]. 
We hypothesize that this operon was highly expressed 
by L intracellularis organisms in the cell endosome 
(Figure 4B). 

We found five genes encoding hypothetical proteins 
among the ten highest expressed by L intracellularis 
(Table 1). From the host gene expression, we identified 



significant activation of Rho family genes which are cru- 
cial for the progression of the Gi phase of the host cell 
cycle and are targeted by bacterial toxins [1]. Taken to- 
gether, this evidence suggests the presence of an 
unrecognized cyclomodulin encoded by those highly 
expressed genes of L intracellularis. Our results form 
the basis for future studies focusing on the specific prop- 
erties and functions of these genes. 

Conclusions 

The present study established a method integrating laser 
capture microdissection and RNA-seq technology to 
characterize both the host and the bacterial gene expres- 
sions in vivo. This methodology was applied to investi- 
gate the proliferative changes induced by an obligate 
intracellular pathogen, Lawsonia intracellularis, in intes- 
tinal crypt cells. The proliferative changes in infected 
enter ocytes were associated with deregulation of the Gi 
phase of the host cell cycle based on the activation of 
transcription, protein synthesis and Rho family genes 
[22,24] . The ability of L, intracellularis in preventing en- 
terocyte differentiation and maturation was proved by 
the consistent down-regulation of apical membrane 
transporters related to nutrient acquisition (Figure 4A). 
The gene encoding copper uptake protein was the only 
apical membrane transporter activated in the Lawsonia- 
infected cells indicating a host defense mechanism by 
inducing copper toxicity [43]. High expression of the 
Cu/Zn superoxide dismutase {sodC) gene by L, intracel- 
lularis suggests that the intracellular bacteria incorpor- 
ate intracytoplasmic copper to cope with oxidative 
stress. In addition to sodC, the intracellular bacteria 
showed high level of expression of a sophisticated oxida- 
tive protection mechanism which involves redox enzymes 
and rubrerythrin-rubredoxin operon (Figure 4B). Finally, 
Rho genes expressed by the host and associated with the 
high expression of bacterial genes encoding hypothetical 
proteins implies a potential role of unrecognized bacterial 
effector proteins that modulate eukaryotic cell cycle 
(cyclomodulins) in the pathogenesis of PE. 

Methods 

Preparation and quantification of L intracellularis 
inoculum 

The porcine pathogenic L. intracellularis strain PHE/ 
MNl-00 was previously isolated and continuously grown 
in cell culture using murine fibroblast-like McCoy cells 
(ATCC CRL 1696) for 10 passages [19]. The bacteria 
were then pelleted, suspended in sucrose-potassium glu- 
tamate (pH 7.0; 0.218 M sucrose, 0.0038 M KH2PO4, 
0.0072 M K2HPO4 and 0.0049 M potassium glutamate) 
solution with 10% fetal bovine serum and stored at 
-80°C until the day of infection. The number of L. 
intracellularis organisms was assessed by direct counting 
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after immunoperoxidase staining [49] and by quantita- 
tive PGR (qPCR), as described elsewhere [50]. 

Experimental infection 

Twelve Duroc-Landrace cross pigs at 3-weeks-of-age 
were divided into two groups: infected and control 
(n = 6/group). The animals were obtained from a herd 
with no history of PE and each treatment group was 
housed in a different pen. Prior to the study, blood and 
fecal samples from all animals were collected and tested 
for L intracellular is -specific antibodies by immunopero- 
xidase monolayer assay (IPMA) and for the presence of 
Lawsonia DNA in the feces in order to confirm their 
negative status [11,49]. The animals were allocated in 
the isolation barns at the College of Veterinary Medicine 
of the University of Minnesota and fed with non- 
medicated nursery feed and water ad libitum. All proce- 
dures were approved by the Institutional Animal Care 
and Use Committee of the University of Minnesota. The 
infected group was orally inoculated with 30 ml of L, 
intracellular is culture at 10^ organisms per ml. The 
non-infected group was orally treated with sodium- 
phosphate-glutamate (SPG) solution. Fecal samples were 
collected every other day and analyzed by qPCR for the 
presence L. intracellular is DNA [11]. 

Tissue collection and preservation 

All animals were euthanized 21 post-infection (PI) ap- 
proaching the peak of the clinical disease, as described 
previously [20]. Immediately after euthanasia, the ileal 
mucosa was washed three times with PBS solution 
containing 2.0 U/(il RNAse inhibitor (Roche Applied 
Science) [51]. A total of six tissue samples from each 
animal were collected, placed into plastic cryomold cas- 
sette (Tissue-Tek® Sakura Finetek) in longitudinal orien- 
tation and embedded in optimal cutting temperature 
(OCT) compound (Tissue-Tek® Sakura Finetek). The 
samples were then placed onto blocks of dry ice until 
the tissues and OCT were frozen to a solid white. All 
samples were stored at -80°C until use. 

Laser capture microdissection 

Five serial 8-(im frozen sections were cut at -20°C using 
RNAse-free blades and mounted on StarFrost® RNAse- 
free slides (Fisher Scientific). The first and fifth 
cryosections (reference slides) were evaluated by immu- 
nohistochemistry (IHC) using the streptavidin method 
with polyclonal antibodies specific for L. intracellular is 
[52]. In these two reference slides the levels of infection 
were assessed based on the amount of positively labeled 
antigen present in the intestinal sections: Grade 0 (-) = 
no positive antigen labeled; Grade 1 (+) = one isolated 
focal area of antigen labeled; Grade 2 (++) = multi-focal 
areas of antigen labeled; Grade 3 (+++) = majority of the 



mucosa has positive antigen labeled; and Grade 4 
(++++) = all of the mucosa has positive antigen labeled 
[17]. The series of slides containing the two flanking 
pairs of reference slides most severely infected (grade 3 
or 4) was selected for LCM from each animal. Similar 
anatomic portions of the ileum were collected from the 
control group and stained by IHC to confirm the nega- 
tive status of the samples. One non-infected series of 
cryosections was selected for LCM from the similar ana- 
tomic portion of those selected in the infected group. 

The LCM procedures were performed in the second, 
third and fourth serially- cut slides using Histogene^^ 
staining (Life Technologies). Briefly, the specimens were 
thawed at room temperature for 30s, fixed in nuclease- 
free 75% ethanol, rehydrated in nuclease-free water for 
30s, stained with HistoGene® solution for 20s and rinsed 
in water for 30s, then dehydrated by sequential 
immersion into 75%, 95% and 100% ethanol for 30s 
each. All microdissection procedures were performed at 
the College of Biological Sciences Imaging Center at the 
University of Minnesota using PixCell® II System (Life 
Technologies). Single-cells were captured from the slides 
using CapSure LCM caps (Life Technologies) with laser 
spot diameter 7.5 (im, power 60 mW and duration 650 
(IS. Approximately, one thousand cells were captured from 
the bottom to the upper part of the intestinal crypts 
(Figure lA and IB) in each of the three stained slides in 
approximately ten minutes, incubated at 42°C for 30 min 
in 50 (il extraction buffer and stored at -80°C. The 
microdissected cells from each of the three serially-cut 
slides were pooled and used for the RNA extraction. 

RNA isolation and amplification 

Total RNA was isolated from microdissected tissues 
using a PicoPure™ Kit (Life Technologies) according to 
the manufacturers instructions. Genomic DNA was 
digested using DNase I (Qiagen). RNA samples were 
then assessed with a NanoDrop ND-8000 Spectropho- 
tometer and Agilent Bioanalyzer 2100. Amplified cDNA 
was prepared from 100 ng of total RNA using the 
Ovation" RNA-Seq System V2 (Nugen"), following the 
manufacturer s instructions. The amplification was initi- 
ated at the 3' end as well as randomly throughout the 
sample which allows the amplification of eukaryotic and 
prokaryotic RNA transcripts. The amplified cDNA 
generated was then purified using MinElute Reaction 
Cleanup Kit (Qiagen) for library preparation. 

Library preparation and sequencing 

The library preparation and sequencing were conducted 
in the core facility of the Biomedical Genomics Center 
at the University of Minnesota. Briefly, 100 ng of the 
amplified cDNA was fragmented, ends repaired, 
enriched by PCR and validated using High Sensitivity 
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Chip on the Agilent2100 Bioanalyzer. Following quantifi- 
cation of the cDNA generated for the library using 
PicoGreen Assay, the samples were clustered and loaded 
on the Illumina® Genome Analyzer GA IIx platform which 
generated on average 22,136,064 paired reads of 100 bp. 
Base calling and quality filtering were performed following 
the manufacturers instructions (lUumina® GA Pipeline). 

Filtering and mapping of sequence reads 

Quality control, trimming and mapping were performed 
in the Galaxy platform [53]. Initially, FastQC tool was 
applied to the raw sequence data followed by FastQ 
Trimmer [54]. As a result, 15 base pairs were trimmed 
from the 5 end and ten from the 3 end of each raw se- 
quence read. Sequences containing 75 bp and a phred 
quality score of more than 20 were used in the gene ex- 
pression analysis [55]. Using Bowtie short read aligner 
with no more than two mismatches, the filtered se- 
quences were mapped onto the pig genome S, scrofa 
10.2 and the L intracellular is isolate PHE/MNl-00 ref- 
erence genome, both obtained from the National Center 
for Biotechnology Information (NCBI) database [56]. 
The number of reads mapped within each annotated 
transcript was calculated in order to estimate the level of 
transcription for each gene. Cufflinks tool was used to 
estimate the relative abundances of the transcript reads 
for each gene [57]. For comparing the levels of gene ex- 
pression between infected and non-infected enterocytes 
and characterizing the abundance of L. intracellular is 
transcripts expressed in the host cytoplasm, the read 
counts were normalized based on the number of reads 
per kilobase of coding sequence per million mapped 
reads (RPKM). 

In order to determine any potential co-infection 
associated with the intracytoplasmic L. intracellular is, 
the sequence reads were screened against 4,279 viral 
sequences available at NCBI database. Additionally, 
MetaPhlAn (Metagenomic Phylogenetic Analysis) was 
performed to identify other bacterial genome sequences 
in the experimental samples [58]. 

Differential gene expression 

The differential gene expression between infected and 
non-infected enterocytes was assessed using the CuffDiff 
tool. This tool defines the total read count for each gene 
by combining the expression data from all the replicates 
in each experimental group and testing for differential 
regulation in expressed transcripts present in at least 
three replicates. RPKM values were expressed in log2 
(RPKM) to allow the statistical comparison. As a result, 
log2-fold change in abundance of each transcript 
was obtained by log2 (RPKM [infected] /RPKM [non- 
infected]). P- values were calculated and adjusted for 
multiple comparisons using false discovery rate (FDR) 



[59]. Significant differential expression was determined 
in genes with FDR-adjusted p-values < 0.05 and fold 
change ± 2. 

Functional gene ontology and pathway analysis 

Biological functions and interactions of the genes differ- 
entially expressed were determined using the database 
for annotation visualization and integrated discovery 
(DAVID v6.7) and the Ingenuity pathway analysis (IPA) 
(Ingenuity® Systems). While the DAVID knowledge data- 
base was used for functional clustering, the mammalian 
knowledge database of the IPA was used for pathway and 
network analyses of the differentially expressed genes. 

Predictive analysis regarding the ontology and motifs 
of gene-encoding hypothetical proteins highly expressed 
by L. intracellular is were performed using Kyoto 
Encyclopedia of Genes and Genomes (KEGG) database. 
Additionally, Phyre^ (Protein Homology/analogY Recog- 
nition Engine) [60] and ConFunc (Protein Function Pre- 
diction Server) [61] software were used to predict their 
biological functions and three-dimensional structure. 

Reverse transcription PGR 

In addition to the assessment of RNA quality based on 
the bacterial and eukaryotic ribosomal RNAs using 
Agilent Bioanalyzer 2100, one-step RT-PCR was applied 
to the RNA extracted from microdissection tissues in 
order to evaluate the quality of bacterial mRNA for 
RNA-seq analysis. The one-step RT-PCR (Qiagen®) was 
used with specific primers (Additional file 3: Table SI) 
targeting three housekeeping genes of L. intracellular is. 

Quantitative RT-PCR (qRT-PCR) was performed for 
validating the expression data identified by RNA-seq 
using a specific set of genes from the host and bacterium 
(Additional file 3: Table SI). RNA samples were synthe- 
sized to first-strand cDNA using Superscript® II RT 
(Invitrogen). Duplicate qRT-PCR reactions from each 
primer probe set were validated by five serial dilutions 
of cDNA on the ABI7900HT instrument (Applied 
Biosystems). After validation, quantitative PCR was 
performed in duplicate using 15ng of cDNA per sample 
with the following conditions: 60°C for 2 min, 95°C for 5 
min and 45 cycles (95°C/10 sec and 60° for 1 min). The 
differential expression for the host gene expression were 
described as relative fold-change between infected and 
non-infected enterocytes based on Ct values (2'^^^^), as 
previously described [62]. The suitability of four porcine 
housekeeping genes (Beta actin, Cyclophilin A, Hypo- 
xanthine phosphoribosyltransferase, Ribosomal protein 
L4) for qRT-PCR analysis was evaluated. 

Based on consistent Cyclophilin A gene expression 
throughout the biological replicates, this gene was used 
to normalize the expression data. For the bacterial gene 
expression, averages of relative transcriptional levels 
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normalized by the SOS ribosomal gene were calculated 
and compared to the RNA-seq expression levels. A lin- 
ear regression model was used to evaluate the correl- 
ation between average RPKM and qRT-PCR data, as 
described previously [63]. 
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